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We constrain the one-parameter (a) class of TeVeS models by testing the theory against both 
rotation curve and strong gravitational lensing data on galactic scales, remaining fully relativistic 
in our formalism. The upshot of our analysis is that - at least in its simplest original form, which 
^\ . is the only one studied in the literature so far - TeVeS is ruled out, in the sense that the models 

cannot consistently fit simultaneously the two sets of data without including a significant dark 
matter component. It is also shown that the details of the underlying cosmological model are not 
relevant for our analysis, which pertains to galactic scales. The choice of the stellar Initial Mass 
Function - which affects the estimates of baryonic mass - is found not to change our conclusions. 
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I. INTRODUCTION 

With the availability of "precision" cosmological data, it has been possible to formulate and test the standard 
(ACDM) cosmological model to unprecedented accuracy, with the paradigm consistently proving itself to be highly 
O !■ successful at fitting observations The model is based on two key components: a framework of classical general 
relativity acting on a homogeneous Friedman-Lemaitre-Robertson- Walker metric with cosmological constant, A > 0, 
and the presence of Cold Dark Matter (CDM). There have been, however, other proposals for models that lie outside 
the ACDM framework. Such models come from those who see the as yet undetected status of dark matter and 
the lack of a fundamental theory for the cosmological constant as impetus for looking into alternative theories. In 
particular, Milgrom Q put forward the Modified Newtonian Dynamics (MOND) to explain the flat rotation curves 
of galaxies without invoking dark matter. Constructed to give flat rotation curves below an acceleration scale, 
ao as 1.2 x 10~ 10 m/s 2 , set by galactic rotation data, MOND modifies the standard Newtonian relation between the 
gravitational potential and the acceleration to f(\a\/ao)a = — V$n, where $n is the Newtonian gravitational field. 
^sO • The function f(x) is positive, smooth, and monotonic; it controls the interpolation between the two gravitational 
regimes above and below the acceleration scale. When f(x) equals unity, the usual Newtonian dynamics holds, while 
I , when f(x) approximately equals its argument, f(x) ~ x, the deep MONDian regime applies. 

MOND has proven itself to be very successful at fitting galactic rotation curves, but it has been less so when 
applied to galaxy clusters. A relativistic field theory counterpart to MOND has also been expounded [3, 4], involving 
the inclusion of a TEnsor, VEctor and Scalar field (TeVeS). Since it has been constructed to reproduce the same 
modification to gravity as MOND in its low acceleration non-relativistic regime, TeVeS has not been tested using 
rotation curves; it has been applied to areas such as the cosmic microwave background [BJ and gravitational lensing, 
| where the relativistic aspects of the theory could be fully exploited in a way unachievable by MOND. Most of the work 
on Te VeS and lensing has remained non-relativistic [6|, |7|, l8| , by simply reducing the theory to that approximately 
replicated through the addition of a scalar-field potential to the Newtonian potential. However, recent work Q 
has conducted a fully relativistic lensing analysis of TeVeS and applied the results to a sample of six lenses from 
the CASTLES survey. In order to ascertain the validity of the TeVeS/MOND claim of explaining observations 
without the inclusion of dark matter, the stellar content of these galaxies, determined by means of a non-parametric 
model-independent method, was compared to the mass predicted from lensing in the context of the TeVeS model. 
A discrepancy here would be indicative of a shortcoming of the theory in its fundamental claim. As the precise 
results of the comparison of these masses depended significantly on the exact form of the MOND f(x) function, and 
its TeVeS equivalent /J,(y), a parametrised range of these functions was considered. The analysis showed [9( that 
the parametrisation commonly used for MOND and TeVeS requires significant quantities of dark matter, with the 
dependence increasing as the parametrisation moved to the functions which best fit rotation curves. These results 
suggest that a harmonisation between fitting both rotation curves and lensing may not be possible in these modified 
gravity theories. 

The purpose of the present work is to fully examine the efficacy of TeVeS, using the two distinct approaches of 
rotation curves and lensing data in order to examine the validity of the implied deficiency found in the theory [||. 
Using the two methods and sets of data, in a fully relativistic way, we independently constrain the degeneracy that 
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exists within the theory regarding the precise form of the n(y) function. In this way we can check whether there is 
any single form of the theory which can fit and describe both rotation curves and lensing data. 

The structure of the article is the following: In Section [TTl we describe the relevant formalism, and solving the 
appropriate classical equations for the metric in TeVeS. We then use this metric to find the modified equations for 
the deflection of light and circular rotation. In Section HTT1 we analyse rotation curve data for a selection of both High 
(HSB) and Low (LSB) Surface Brightness galaxies. In Section HVl we analyse gravitational lensing for the constraints 
imposed by the rotation curves. We also find the parameter which fits best the lensing data, and examine how this 
fits with rotation curves. We round up with an analysis of the results, and present our conclusions and outlook in 
Section El 

II. THE TEVES MODEL 

In this section we shall review the relevant formalism, leading to the TeVeS dynamical equations, which are used 
when computing the deflection angle and other quantities to be used in our fitting procedure. We shall be brief in our 
analysis. For details we refer the reader to reviews and relevant previous work [a. M. IToL lllj] . TeVeS Q is a bi- metric 
model in which matter and radiation do not feel the Einstein metric, g a p, that appears in the canonical kinetic term 
of the (effective) action, but a modified "physical" metric, g a p, related to the Einstein metric by 

g aP = e- 2 *g a(3 - U a Up(e 2 * - e" 2 ^) , (2.1) 

where U^, <p denote the TeVeS vector and scalar field, respectively. The TeVeS action reads 



d 4 x 



' (R - 2A) - \{a 2 (g^ - £W) <f> ia <P, + ^Gl^^FikGcj 2 )} 
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{-g) 1/2 +C{g, v ,r,f^,...){-~gfl 2 , (2.2) 



where k, K are the coupling constants for the scalar, vector field, respectively; £ is a free scale length related to ao 
(c.f below); a is an additional non-dynamical scalar field; T ^ = Un jt/ — U v „; A is a Lagrange multiplier implementing 
the constraint g a ^U a Up = —1, which is completely fixed by variation of the action; the function F(kGcr) is chosen 
to give the correct non-relativistic MONDian limit, with G related to the Newtonian gravitational constant, Gn, by 
G = Gn/[1 + K/2 + fc/(47r) — 2<p c ], where <j> c is the present day cosmological value of the scalar field. Covariant 
derivatives denoted by | are taken with respect to g^v and indices are raised/lowered using the metric g^. 

The modified equations of motion can be calculated from the Lagrangian. For the modified Einstein equation we 
have 0,[lil 



Gap + .9c/3A = 8itG 



T a p + (1 - e-^)U»T Ka U p) + r afj + Q a p , (2.3) 



where 



t q/3 = a 2 [(j)^ a (j)^ ~ -g^&^^gap - U^(j) tfi (U^ t ^ - \u v <\>, v g a p)\ - ^ T F(kG(7 2 )g a p , 

Q af3 = K(g^U M U M - \g aT g^U [a ^U lTM g a p) - \U a U . (2.4) 

For the vector field we obtain 

8ttG(1 - e-^g^U^f^ = KU [a ' J3] , p + XU a + 8ttG(7 2 C/ /3 f <f>,pg ai '<j> n , (2.5) 
and similarly for the scalar field, namely 

Hv)(9 a0 - U a U^ a y = kG{g a ? + (1 + e-^)U a U?]f a p , (2.6) 

with /i(y) defined by: 

H(y) = kGa 2 , 

1 2 dFV) 

y = -mF(m) - ^ , 

y = kl 2 (g^ - U^U v )<t>,^ v . (2.7) 
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Motivated by the homogeneity and isotropy of the Universe, observed to date, we assume a spherically symmetric 
metric 

g a pdx a dx p = -e u dt 2 + e c {dr 2 + r 2 d6 2 + r 2 sin 2 9d<p 2 ) , (2.8) 
where v and £ are both functions of r. The physical metric can be written likewise, namely 

g a pdx a dx f} = -e^dt 2 + e<(dr 2 + r 2 d6 2 + r 2 sin 2 6»d^ 2 ) , (2.9) 

with the quantities v and £ related to v and £ by 

v = v + 2<j) ■ C = C-20. (2.10) 

Isotropy makes the scalar field depending only on r, namely cf> — <fi(r). We approximate matter as an ideal pressure- 
less fluid, T a p = pu a up. By assuming that the time-like vector field has only one non-zero temporal component, the 
normalisation condition imposed by the Lagrange multiplier restricts the vector field to be 

U a = {e- v/2 , 0,0,0) . (2.11) 

Considering the quasi-static case, we can take the four- velocity of the fluid, u ai to be collinear with U a , and then 
normalise it with respect to the physical metric, g a p, so that u a = e^U a , leading to 

f a() = pe 2 *U a U . (2.12) 
Thus, the scalar field equation, Eq. (|2.6j) . along with the isotropy constraint, leads to 

rVe^^(y)l ' = hGe~ 2 *p , (2.13) 



where a prime denotes derivative with respect to r. Upon integration, we obtain 

_ kGm a (< r) „_ (v+c)/2 



4irr 2 p(y) ^ ^ 



where a scalar mass has been defined as 



m s (< r) — 4ir / pe^ + ^ 2< V 2 dr . 
Jo 

As shown in Ref. [|| , the scalar mass can, to a good approximation, be equivalent to the "proper" mass contained in 
the same volume. Moreover, the Lagrange multiplier appearing in the vector field, Eq. (|2.1ip , can be totally determined 
by the vector equation, Eq. (|2.5p . namely 



A = 8^G(e" 2 * - e 2 *)p - Ke~< (y + ^ + 7) 



(2.15) 



Solving the modified Einstein equations, Eq. (|2.3p . for v and £ we find for the (tt) and (99) components of the stress 
tensor 

kGm s (< r) 2 p 2 (y) v 

" ~ 327T 2 /x(y) r 4 + 4GZ 2 fc 2 W ' 



2 ft -(C+«) ,,2 



Tee = 



i/-n2i.2 x W e ' 



kGm s (< r) 2 e 
327r 2 ^(y) r 2 AGl 2 k 2 

(2.16) 

0« = Ke»-t (^f- + ^ + + ^\ - 8nGp(e- 2 * - e 2 *)e» , 

®ee = ^(rv'f . (2.17) 
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Moreover, the Einstein tensor components are: 
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We thus arrive at the following system of differential equations: 



(2.18) 
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We analyse two distinct systems: lensing galaxies and galactic rotation curves. For the former, we model the system 
with the spherically symmetric mass profile [l2T |: 



M{< f) 



Mr 2 



(f + r h ) 2 



(2.20) 



where f 



oC/2 



r is the Schwarzschild radial coordinate, is the core radius scale, related to the projected two- 



dimensional half mass radius, R e , by R c = 1.8153rh [13|, and M is the total mass of the galaxy. The mass, Eq. (|2.20[) . 
specifies the density function p, and, as we mentioned previously, it is assumed approximately equal to the scalar 
mass m s (< r). 

For the galactic rotation curves, we use a different spherically symmetric mass profile, as outlined in Ref. [l4| . given 
in Schwarzschild radial coordinates by 



M(< f) = M 



3£f 



(2.21) 



where (3=1 for HSB galaxies and (3 = 2 for LSB galaxies; r c is the core radius and M is the total mass of the galaxy. 
We remark at this point that, when converting to the physical coordinates, a factor e^ +2 ^ appears as the coefficient 
of the cosmological constant A. 

In order solve the equations analytically, we have used the following approximation Q: 



e 20 ^ e 20 c 



kGm s (< r) k 2 G 2 m 2 (< r) 



2irr 



8jT 2 r 2 



Q{r- 



(2.22) 



The precise form of the modification to gravity given in TeVeS is largely controlled by the fi(y) function, however, 
since the theory is not motivated from a microscopic theory, there exists a large amount of freedom in choosing the 
form of this function. In Ref. [|[ a toy-model function was suggested, namely 



y 



3 ^ 2 (^-2) s 
4 1- n 



(2.23) 



However, it was noted in Refs. [15| and [l6( that when this function is converted to its MONDian equivalent, it fits 
rotation curve data worse than the standard MONDian "simple" form 



/(*) = 



In order to improve on this point, the authors of Ref. [l7| suggested the following function 



fix) 



2x 



1 + x + - x) 2 + Ax 



(2.24) 



(2.25) 



which fits better the rotation curve data in the MONDian framework. 
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In addition, they provided parametrised functions f(x) and n{y) which interpolate between the toy function, 
Eq. (|2.23[) , and the "simple" function, Eq. (|2.24p , in both MOND and TeVeS frameworks, over the parameter range 
< a < 1: 

2x 

fix) = = ; u(y) = Y 3 p= ; < a < 1 . (2.26) 

W l + (2-a)x + ^/(l-ax) 2 +Ax ^' 1-2^^' " 1 ' 

Note that, for a — 0, the function n(y) does not coincide with the the toy function, Eq. (|2 . 23[) . but rather approximates 
the function in the high and intermediate gravity regimes. However, this approximation is valid for our present analysis. 

In fact, for our purposes, we find it convenient to use an explicit parametrisation of /i in terms of the scalar mass 
m s , which is approximated by the mass profiles given directly by the data. To this end, we use Eq. (|2.14p and the 
definition of y, Eq. (|2.7p . as well as the scalar equation of motion, Eq. (|2.6p . to write the parametric function fi(y) in 
Eq. (|2T2"rJj) as: 

2naj 2tt fa 2 ] 2 k 2 j \ 1/2 , , , / fc 3 / 2 IG\ m s (< r)e^l 2 -^ , n ^ 

The corresponding F function is then given by (c.f. Eq. (|2.7[) ~): 

6fc 3 



(47ra) 3 /i 2 



j / Ana/A ^\ " \ Aira/i 



\ k J 1+ ^LBE 

k 



(2.28) 



When a = 0, the function F((a) becomes singular. Since the a = case is supposed to give a very close approximation 
of the toy model function, we use for a — the explicit expression for F(n) given in Ref. [3J|, which does not become 
singular. With all the various components given above, we are able to numerically solve the system of differential 
equations to find the functions £, v, which specify the metric for our choice of the /i(y) function. 



III. FITTING ROTATION CURVES IN TEVES 



^ = ^=(-^1 (3.1) 



In our metric system, Eq. (|2.9[) . the equation for circular orbital velocity in TeVeS reads 

1/2 

d< ~ \2 + rC 



In order to perform a fitting analysis to rotation curve data there are two free fitting parameters: r c , which is the 
core radius of the galaxy, and a, the parameter appearing in the TeVeS fi(y) function. For a particular choice of 
these parameters, we can predict the expected rotational velocity at certain radii and this is compared against data. 
Moreover, the constants in the TeVeS action, Eq. I|2.2p . are taken to be 



= 0.01 ; if = 0.01; £ = Vkh / {AnEao) ; 
cf> c = 0.001 ; 3 = 1 + K/2 - 2(j) c . (3.2) 

The values of K and k are constrained 3] from solar system tests on gravity to be < 0.1, and by rotation curves to be 
> 0.001. The scale I is related to the MONDian acceleration scale, a and b. The latter quantity is found by taking 
the limit of the function y(/i) when /i << 1, which then takes the form y([i) ~ so for the class of fi functions 
considered here, we set b — 3. Finally we note that, for cj) c , the present day value of scalar field at cosmological scales, 
there are no tight constraints on its exact value, with an approximate upper bound coming from cosmological data. 

With these constants so defined, we are able to numerically solve the equations and perform a fitting analysis for 
the galactic rotation curves. A standard likelihood analysis is used to determine the best fit parameters. Using x 2 - 
statistics, the likelihood probability distribution function for the value of a using data from the i th galaxy is defined 
as 

Li(a) = e-ilXiM-x'n.iJ . (3.3) 
The combined likelihood function for all the galaxies is therefore defined as: 



Ltot(a) = Lx{a) x L 2 (a) x L 3 (a) x • • ■ , 



(3.4) 
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M K 


M gas 
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;amag 






xlO 8 M 




xlO lo M 


NGC 3972 


0.97 


0.43 


1.60 


-22.08 


18.07 


57.15 


80.17 


3.77 


NGC 4085 


0.99 


0.46 


1.70 


-22.27 


15.89 


72.09 


100.93 


2.75 


UGC 6399 


0.88 


0.33 


1.56 


-20.33 


11.43 


11.30 


15.85 


1.46 


UGC 6917 


0.86 


0.35 


1.26 


-21.10 


28.52 


27.61 


38.37 


3.10 


UGC 6923 


0.83 


0.53 


1.16 


-20.36 


11.65 


13.96 


19.41 


0.84 


UGC 7089 


0.85 


0.30 


1.04 


-20.30 


18.51 


13.77 


19.14 


1.37 



TABLE I: Photometric and baryonic mass data for the galaxies in our sample. The photometric data are corrected for 
internal extinction and are taken from from Table 2 of Ref. M gas is the gas mass from 21 cm fluxes corrected to include 

Helium, M st is the calculated mass of the stellar content (given for both choices of the Initial Mass Function). Mn gives the 
Newtonian estimate for the total mass, by using the maximum observed radial position and the maximum rotation velocity: 
Mn = «max-R/G. 

where the subscript on L denotes the galaxy. The peak of the total likelihood represents the most likely value of a, 
and the cumulative distribution is used to determine the 90% confidence level (i.e. comparing the 5- and 95-percentile 
of the distribution). The constraint on a from this analysis is then used to determine lensing masses, as shown below. 

In order to find the minimum value of x 2 by numerically varying the free parameters, an adaptive grid method was 
implemented to give a value of \ 2 accurate up to 2 decimal places. The sample of galaxies we chose to examine were 
taken from Refs. [HI, [lj| which explored the Tully-Fisher relation in galaxies belonging to the Ursa Major cluster. 
Only those rotation curves were examined for which there exist at least five data points and the inner part of the 
velocity distribution is available, i.e. we excluded those galaxies for which the data only recorded the velocity in the 
flat part of the rotation curve and missed the inner sections. This gave us a sample of 6 galaxies, 2 LSB (NGC 3972, 
NGC 4085) and 4 HSB (UGC 6399, UGC 6917, UGC 6923, UGC 7089). The photometric details of the galaxies are 
given in Table I. 

The photometric details given in Table I are used to calculate the total baryonic content of the galaxies, fixing the 
mass parameter for TeVeS as M gas + M st . For the gas mass, we assume an infinitely thin disk and we use the observed 
21 cm line fluxes (listed in Ref. [l8j), using the standard translation between 21 cm flux and neutral hydrogen mass 
(see, e.g. Ref. [l4j). A correction factor of 4/3 is necessary to take proper account of the presence of Helium in the 
gas phase. 

For the stellar content of the galaxies it is important to choose realistic populations. Population synthesis models 
(e.g. Ref. [20]) combine our knowledge of stellar evolution with libraries of stellar spectra to obtain simple stellar 
populations (SSPs), which correspond to a family of stars all formed at the same time, with the same metallicity. 
Although SSPs can be used to describe the spectral energy distribution of globular clusters - which form over very 
short times compared to stellar evolution timescales - a galaxy has a more complex distribution of stellar ages. We 
follow here the same approach as in Ref. [j~3l ] , and run a grid of models with an exponentially decaying star formation 
history. These models replace the single age of a SSP by a linear superposition of SSPs according to an exponentially 
decaying function of time. Hence, a model is described by the following three parameters: (i) the epoch when the 
galaxy is born (chosen between redshifts zp = 1 and 10), (ii) the timescale of the exponential (between logr(Gyr)= — 1 
and +1) and (iii) the metallicity (between l/30th and twice the solar abundance). Our "basis" SSP models are taken 
from Ref. [21]). We run a grid of 32 x 32 x 32 models with a uniform distribution of the three parameters that describe 
the star formation history. For each choice, the models give a spectral energy distribution which is folded with the 
response function of the standard filters B, R, I and K, which extend over a wide spectral range, from optical blue 
light (A ~ 0.4/im) to near infrared (A ~ 2.2/xm). Our photometric system is referenced with respect to Vega. The 
resulting colours B — R, R — I and I — K are compared with the observations, defining a likelihood which is used to 
determine the stellar masses by comparing the best fit absolute magnitude in the near infrared band (K, A ~ 2.2/im). 
At longer wavelengths, stellar mass estimates are more robust, since they minimise the attenuation from dust and the 
flux mostly originates from intermediate/low mass stars, which dominate the stellar mass budget. 

In addition to age and metallicity, the stellar mass of a population depends sensitively on the mass distribution 
of stars at birth, i.e. the Initial Mass Function (IMF). It is generally assumed that the IMF is universal, although 
its shape at the low mass end is uncertain given how little low-mass stars contribute to the total luminosity of a 
population. In this paper we consider both the Chabrier IMF [12] and the simple power law of Salpeter [23|. Both 
IMFs roughly agree at the massive end - where observations can set significant constraints - however, at the low mass 
end, the Salpeter IMF simply extrapolates the distribution as a power law. Recent dynamical and lensing studies 
(see, e.g. Ref. [25|,[26|]) suggest that the Salpeter IMF gives unrealistically high stellar masses, and functions such as 
the Chabrier IMF - which replace the power law by a log-normal distribution at low masses- are favoured. However, 
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UGC 6917 


9 


13.84 


1.12 


20.91 


9.56 


1.23 


17.59 


UGC 6923 


4 


11.66 


1.31 


4.8 


12.14 


1.30 


4.87 


UGC 7089 


10 


8.53 


0.58 


2.71 


7.12 


0.6 


2.26 



TABLE II: Details of the TeVeS parameter fits for the rotation curves, given for both choices of the Initial Mass Function: 
Salpeter and Chabrier . Note, that the degrees of freedom are denoted by d.o.f. When the best fit parameter is a — * oo, it 
represents cases where the best fit is given by GR (TeVeS tends to GR as a goes to infinity) . In this paper we also discuss the 
a — case (which gives the lowest Dark Matter fractions for strong lensing in TeVeS. For comparison we include the reduced 
X 2 for this case. 



given that all those analyses were based on a "standard" General Relativity (GR) / Newtonian theory, we present in 
this paper both choices of IMF (see, also Ref. (27|). 

The total baryonic content is the sum of the stellar and the gas mass. Within the TeVeS framework, we assume 
that the baryonic masses are the total masses of galaxies and we combine them with the observed rotation curves 
to constrain the a parameter that characterises TeVeS. The values of the likelihood probability distribution function 
fits are plotted against the value of the a parameter. The results are given in Fig. [1] as a two-dimensional (2D) map 
with a vs. core radius. The contours are shown at the 75%, 90% and 95% confidence level for the Chabrier (red) and 
the Salpeter (blue) IMF. We also show the marginalised ID likelihood function with respect to a in Fig. [2) with the 
results for the Chabrier and the Salpeter IMF given as a solid and dashed line, respectively. The values for the lowest 
X 2 and the corresponding values of a are given in Table II, where it can be seen that each of the galaxies analysed 
here gives a reasonable fit to the data once the x 2 is minimised. The results show that the best fit values for the a 
parameter are in the range ~ (5 — 10). 

Figure [3] shows the result for the combined probability distribution function over the six galaxies. By normalising 
the total probability distribution function and setting confidence limits at the 95% level, the uncertainty for the best 
fit value is calculated: a = 11.56t4 2 8 Q 7 f° r ^ ne Salpeter IMF and a = 8.54^3° 32 for the Chabrier IMF. The constraints 
on the a parameter are then used in an analysis of strong gravitational lensing on galaxy scales in TeVeS, to assess 
the compatibility of these two different ways to constrain mass over similar length-scales. 



IV. GRAVITATIONAL LENSING IN TEVES 



In this section we perform a gravitational lensing analysis of the TeVeS models. Given the form of our metric, 
Eq. (|2.9p , we can derive the equation for the deflection of light in the physical metric. It reads @ 



A0 



1-1/2 



b 2 



Ar — 7r , 



(4.1) 



where b is the distance of closest approach for the incoming light ray and it is related to r , the impact parameter 
through 



^{ra)-Kro) r 2 



(4.2) 



Once more, these equations can be numerically solved using the same definitions of the TeVeS constants given earlier. 
In order to perform the lensing analysis in TeVeS, we use the result, Eq. (|4.1[) . for the deflection angle in the lensing 
equation (see e.g. [28|) 



a{9, M, b) 



As 



(4.3) 



which relates the actual position of the background source (3 to the observed position of the images, given by 0. For a 
given cosmological model, the angular diameter distance from the lens to the source, D\ s , and from the observer to the 
source, D Sl are both taken from the redshifts. The lensing equation is applied independently to the multiple images of 



8 




Core Radius (kpc) 



FIG. 1: Likelihood maps of the best fit for the rotation curves of disk galaxies, treating the core radius (horizontal axis) and the 
TeVeS a parameter (vertical axis) as free parameters. The red (blue) contours correspond to the 75%, 90% and 95% confidence 
levels of the Chabrier (Salpeter) Initial Mass Function, respectively. 



the background source in order to find the actual solution of the lcnsing model, which gives the actual position of the 
source (/3) and the mass of the lens. Given that these lenses arc early-type galaxies, a Hernquist profile [12] is adopted. 
This profile has been chosen because its projection gives the characteristic i? 1 / 4 surface brightness profile of this type 
of galaxies (see Ref. [HJ for details). We adopt a concordance cosmology, namely (f2 m ,^A,^k) = (0.3,0.7,0). As 
shown in Ref. @, deviations from this cosmological picture lead to very little difference in the analysis. 

Figure 2] illustrates the methodology we use to determine the lensing masses, for lens HE1104-1805 as an example. 
The lensing equation, Eq. (|4.3[) is represented by the pairs of curved lines that intersect at the true value of the lens 
position and lensing mass. Three pairs of lines are shown for the best fit case of TeVeS (solid black); standard GR 
(grey solid) and the a — case (dashed), which would be the best choice for a lensing based study, although as shown 
above, it is significantly ruled out by the rotation curve data. The figure shows the result both for a Chabrier (left) 
and a Salpeter IMF, with the stellar masses given by the big solid dot and the error bars. 

To calculate a possible dark matter component, a sample of double lensing systems from the CASTLES database [HJ 
is analysed and the mass of the lensing galaxy in GR and TeVeS is calculated, using the same fully relativistic method 
outlined in Ref. @ . By comparing the mass from lensing to the stellar mass content calculated from a comparison of 
photometry and stellar population synthesis using a Chabrier Initial Mass Function (IMF), as in Ref. [13], the mass 
deficit which belonged to the "dark" sector is found. 

From the rotation curve analysis, the favoured values of a are 11.56 and 8.54 for Salpeter and Chabrier IMFs, 
respectively, and these values are examined in the context of lensing. The results from the lensing analysis are given 
in Table III. 

From Table III one can clearly conclude that when a = 11.56 and a = 8.54, there is a considerable need for dark 
matter, comparable to that required in the pure general relativistic case. We find that on average over the six galaxies, 
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0.81 
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+0.05 
0.06 
+2.10 
2.36 
0.98 
0.41 
0.29 
0.35 



30.0 
47.2 
53.2 
49.4 
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1.8 

2.1 
+0.9 
1.2 
+1.2 
1.5 
+3.2 
1.5 
1.1 
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14.40 
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+0.36 
0.45 
+0.28 
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58.6 
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42.1 


21.0 


47.1 


51.6 


22.3 


8.7 


36.8 


14.9 
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TABLE III: M„ is the stellar content of the lenses taken from [251 ] . The other columns show the masses of the lensing galaxies 
for GR and TeVeS when a = 11.56, a — 8.54 and a = 0. It is also given is the percentage of dark matter needed to account 
for the lensing observed given the luminous content of the galaxies. All masses are given in units of 1O 1O M . The errors for 
the a = 11.56 and a — 8.54 are derived from the errors on the value of a and are given at the 95% confidence level. 




FIG. 2: Marginalised likelihood distributions with respect to the TeVeS a parameter for all six disk galaxies. The solid (dashed) 
lines correspond to the Chabrier (Salpeter) Initial Mass Function, respectively. 



in the Salpeter case, TeVeS only predicts 3.4% less dark matter than the general relativistic case: GR requires 47.5% 
dark matter as opposed to TeVeS which requires 44.1%. For the Chabrier case TeVeS gives on average 59.6% dark 
matter as opposed to GR which gives 62.6%, a reduction of only 5.0%. Using the errors on the value of a, the errors 
on the mass are calculated at the 95% confidence level. The error bars are quite tight, only 3.3% on average in the 
Salpeter case and 2.2% for Chabrier, showing just how strongly the analysis rules out this class of fj,(y) functions. 
This result is in agreement with the work detailed in ref. , where it was shown that even for a value of a as low as 
there is a not-insignificant need for dark matter. However the analysis presented in Ref. 1 shows that only values 
of a greater than or equal to 1 can be definitively ruled out by the lensing data due to the error bars on the stellar 
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FIG. 3: Total marginalised likelihood distributions with respect to the TeVeS a parameter combining the data for all six disk 
galaxies. The result for the Chabrier (Salpeter) IMF is shown as a solid (dashed) line. 



masses. The a = case is not ruled out as the lensing masses are within the error bars of the stellar masses, although 
they are consistently just within the upper error bar which suggests that with higher accuracy stellar mass estimates, 
even this case could be ruled out using just a lensing analysis. In the absence of more accurate estimates, we can 
test the a = with rotation curves. Table II comparison of the x 2 for rotation curves for different values of a. As 
it can be clearly seen, for three of our sample of six galaxies (N3972, N4085 and U6917), the a = gives a value of 
X 2 which is nearly an order of magnitude larger than the best fit a case for both Salpeter and Chabrier IMFs, and 
is sufficiently high so as to be considered a poor fit. This implies that the a = case is incompatible with rotation 
curves. This result was also suggested in Refs. [TBI [l6j where the authors showed, using a non-relativistic photometric 
fitting approach to rotation curves as opposed to the relativistic parametric fitting used here, that galaxy NGC3198 
could not fit the data when a = but only when a = 1, a case already ruled out by the lensing analysis of Ref. [s|. 

From this we can infer that there is no value of a which can satisfactorily fit both rotation curves and gravitational 
lensing. In conclusion, at the very least an entirely different form of the /x(y) function, and its related MONDian 
equivalent f(x), needs to be found if the modifications to gravity are to remain universal in applicability with only 
the acceleration scale being free to be fitted. 



V. CONCLUSIONS AND OUTLOOK 



In this paper we have made an attempt to constrain TeVeS models by using both gravitational lensing and rotation 
curve data. In particular, we have analysed the one-parameter models that have been fitted against the rotation 
curves, and we have found that there is no value of the parameter a that fits both the rotation curves and the 
gravitational lensing data of galaxies. The baryonic mass of the galaxies is calculated using photometric data and 
is assumed to account for the total mass budget of the system within the TeVeS paradigm. A standard likelihood 
method gives a = 11.56±4 2 8 o 7 for the Salpeter IMF, and a = 8.54±gf for the Chabrier IMF, at a 95% confidence 
level. Consequently, we estimate the mass content of five strong gravitational lenses from the CASTLES survey 
and compare their lensing masses to the corresponding stellar content, calculated from photometry. On taking into 
account the constraint from rotation curves, we find that the lensing mass within the TeVeS formalism still shows an 
excess around 50% over the baryonic content. The only successful parameter value from lensing (a = 0) is shown to 
be incompatible with rotation curves. 

For the fits we used a particular set of galaxies, for which information for the inner part of the velocity distribution 
is available. The upshot of our analysis is that, at least in its simplest original form, which by the way is the only 
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FIG. 4: Constraining the lensing mass of galaxy HE1104-1805 (z= 0.73). The lines come in pairs, and represent the constraint 
on true source position (/3, vertical axis) and lensing mass (horizontal axis) from each of the two observed positions of the 
images of the background source. The intersection of each pair of lines give the lensing mass. The standard GR solution is 
shown as a grey solid line, the best fit is shown as a black solid line and the a — case discussed in the text is given as a 
dashed black line. Each panel corresponds to a different choice of Initial Mass Function, as labelled. The stellar mass is given 
as a large black dot along with error bars, from Ref. [l3| . 



one studied in the literature so far, TeVeS is ruled out. One may think that multi-parametric extensions of TeVeS 
models may be successful. However, such models have not been proposed so far. Lacking a fundamental microscopic 
derivation of the /i-functions, the use of such complicated models becomes unnatural and probably of no physical 
relevance. 

Several features of the TeVeS model cannot be excluded at present, and in fact such features may be desirable 
and essential for fitting extended versions of TeVeS which include dark matter components. In particular, the bi- 
metric nature of TeVeS is an important ingredient which may indeed characterise more fundamental models, such as 
string- inspired foam ones [!§]. Moreover, the vector- like cosmological instabilities implied by the vector time- like field 
may play an important role in structure formation, as emphasised in Ref. [3fj| . reconciling the observed baryonic 
spectrum with theoretical predictions. It is worthy of mentioning at this stage that in the context of the string 
models of Ref. [29| . the bi-metric feature of TeVeS is realised because of the presence of non-trivial dilaton fields, 
while the vector instabilities arise because of physical interactions of open string matter with space-time defects that 
such models include, as a result of the recoil of the latter. The models have a non-trivial dark energy component built 
in, as a result of the Born-Infcld-type dynamics of the vector field, that is characteristic in open string theories on 
brane worlds considered in that work. In addition, the string spectrum contains dark matter particles, in particular 
supersymmctric ones, and in this sense the TeVeS-like features co-exist with extra components of dark matter and 
dark energy in the model, fitting perfectly all currently available data, including cosmological ones. We should stress, 
however, that in such string models, the fundamental physics is entirely different, and even if one has TeVeS-like 
features, such as bi-metric models and vector-like instabilities, these are features that pre-existed the specific TeVeS 
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models of Ref. [3| and their origin is traced back to fundamental structures in string/brane framework. 

Also, for reasons explained in Ref. (29j . an important role is played by neutrino matter in such models, which 
provide a significant component of dark matter, in addition to that offered by the supersymmetric matter. In this 
latter respect, we also mention other phenomenological works within the context of the TeVeS model or extensions 
thereof [5|, l31(, in a conventional field theoretic framework, claiming a prominent role of neutrinos as dark matter 
components in TeVeS-like models, fitting cosmological data. These are interesting avenues for research, which we 
intend and hope to be able to pursue in the near future. 
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